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ABSTRACT 


An Acousto-Optic Spectrometer (AOS) used as an IF spectrometer to a 
heterodyne receiver Is modeled as a total power multi-channel Integrating 
receiver. Systematic noise contributions common to all total power, time 
integrating receivers as well as noise terms unique to the use of optical 
elements and photo-detectors In an AOS are Identified and discussed. In 
addition, degradation of signal-to-noise ratio of an unbalanced Dicke receiver 
compared to a balanced Dicke receiver is found to be due to gain calibration 
processing and is not an instrumental effect. 
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FOREWORD 


The advantages of using acousto-optic techniques as IF spectrometers come 
from the simplicity of a small number of discrete components necessary to 
construct a system. This simplicity gives rise to compactness, low power 
requi rements , and low weight which Is even more important as heterodyne 
receivers are pushed Into frequency ranges where remote platforms are necessary 
to minimize signal attenuation from the earth's atmosphere. In addition, the 
use of waveguide and integrated-optic components offers an additional degree of 
miniaturization which can further simplify alignment and ruggedization problems 
of remote IF spectrometers. 

In this report we will enumerate possible sources of noise which will be 
associated with the use of acoufito-optic techniques as total power radiometers. 
Some noise sources are common to all multi-channel power radiometers while other 
noise sources are unique to the use of optical components and photodetectors as 
essential elements of the radiometric system. We will see that additional 
signal to noise improvement can be achieved by digital processing over analog 
integration on the photodetectors. However, the use of analog to digital 
converters, ADCs, requires careful consideration of the linearity of the data 
conversion to insure that one can obtain the necessary sensitivity to small 
signals required for astronomical observations. In addition, we derive a 
degradation in the signal to noise ratio that occurs when instrumental 
calibration is taken Into account in an unbalanced Dicke-switched system. 

Noise source measurements in both bulk and integrated optics AOS's will be 
presented in a future report. Tms study will form the basis of a systematic 
comparison between actual devices. 


I. INTRODUCTION 


Whether we use bulk or Integrated-optics techniques, an acousto-optic 
spectrometer (AOS) Is composed of four main components; laser light source, 
collimating and focusing optics, acoustic transducer, and detector array. The 
principle of operation remains the same for both types of system. A conceptual 
diagram of an AOS is shown in Fig. 1. An ultrasonic delay line is used to 
convert wi de-band RF input to a proportional spatial pattern by means of a 
traveling pressure wave in a crystal. The spatial variation created in the 
index of refraction of the crystal is used to modulate coherent light, and the 
diffracted spectral components are separated by a lens. The Fourier transform 
of the RF input signal is distributed on the output plane of the lens, and the 
light intensity is detected by an array of photodiodes, yielding the 
instantaneous power spectrum of the RF input signal. 

Specifically, a given RF signal, g(t), can be decomposed by the Fourier 
transform relationship into plane wave components (cf. [I])- 


g(t) 



G(f)e' i2irtt df 


U) 


where G(f)df can be interpreted as the amplitude of a plane wave of frequency f 
launched by an ultrasonic transducer into a crystal , or onto a substrate surface 
by a SAW transducer as in an integrated-optics AOS. 

The interaction of this wave train (equation 1) with the optical beam 
results in a superposition of plane-wave deflected beams provided that the Bragg 
condition is satisfied for each component, i.e. sinag « a/2a where x is the 
wavelength of light, and a is the acoustic wavelength. 


Often it is necessary to design a set of phased transducers to achieve 
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wideband RF response since It Is impossible for a single transducer to meet the 
Bragg condition for all Input frequencies. The propagating acoustic wave can be 
described by replacing t with t Q 


9(t) 



(2irf/v)z 

df 


( 2 ) 


where t Q is the reference time (given by time of excitation), v is the acoustic 
velocity, and z is the coordinate aligned with the direction of the acoustic 
wave vector. 

The deflected optical beam generated in a Bragg cell modulator is 
proportional to the product of the incident beam wavefront E^(z}, the acoustic 
wavefront g(t), and the efficiency of the interaction n(0- A fraction of the 
light acquires a new direction of propagation as a result of the Interaction, 
with the angle if deflection given by a»xf/v (see Fig. 2). 

The deflected beam E d (z) emerging from the acousto-optic interaction is 
given by 


F-d(z) 


J E j (2)G(f)n(f)e 1 ( 2 l,f ' /v)z dr' 


(3) 


In coherent illumination, a lens of focal length F will give at its back 
focal plane a field distribution which is the Fourier transform of the field 
distribution in its front focal plane; i.e. 


E 2 O0 » c 



e i(2ir/xF)yz d5 


( 4 ) 


where C is a constant, and y is the transverse coordinate axis, in the back focal 
plane. If E d (z) is not in the front focal plane, the output distribution 

is not changed, but a constant phase factor e* x (^ is added to equation 
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FIGURE 2. 


Optical wavefront propagation through the acoustooptic 
interaction region and the transform lens. 
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(4) that is not detectable in power measurements. 

Combining equations (3) and (4), and substituting with 2 cos a to account 
for the new propagation direction gives the field distribution at the detector 
plane 


My)> 



1 2ff (zcos$/v)f i(2*/xf)yz 
G(f)n(f)w(z)e e dfdz 


(5) 


where w(2) is a weighting function determined by the optical aperture size, 
Incident optical wavefront distribution, and system aberration. For example, a 
Gaussian laser beam and an aberration-free system will yield a truncated 
Gaussian form for w(z). 

Carrying out the integration over z in equation (5), and denoting W as the 

» 

Fourier transform of w(2) gives 

E 2 (y)-C J G ( r 'M f ) w ( y + F*cos«/v)df (6) 

with y + (Fxcosa/v) the argument for W. 

We see from equation (6) that the output intensity distribution |E 2 (y)|^ Is 
positioned linearly in proportion to the frequency of the modulating signal and 
its intensity is proportional to the strength of the spectral component at 
frequency f multiplied by the efficiency of the Bragg modulator. A well 
designed optical spectrum analyzer should have frequency as well as intensity 
linearity. 

The intensity of the diffracted light in the Bragg region is given b> [2] 


I, - n Sine I, 2 (n+ Uk L ) 2 ) 1/2 

r i 2 (7) 
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where 1^, and l Q represent the diffracted and Incident light Intensity, ak the 
momentum mismatch between the Incident light and acoustic propagation vectors, 
and L Is the Interaction length of the acoustic wave column 


r, - 




( 8 ) 


where p Is the photoelastic coefficient, p Is the mass density, v Is the 
acoustic velocity, H is the height of the acoustic beam, P a is the acoustic 
power, and the quantity n 6 p 2 /pv 3 Is the acoustic figure of merit Mg. In the 
small argument limit equation (7) yields 


*1 ■ n ' . (9) 

»o 

so that n gives directly the diffraction efficiency of the device. 

For an acoustic modulator with RF bandwidth At', the total angular 
deflection as * af x/v, the number of resolvable angles N Is determined by the 
total angular deflection divided by the resolvable angular spread, X/L, of the 
light beam* 


N - - * Af . ,_L » Af L » tAf (10) 

X/L v x v 

where L/v * T is the aperture, i.e. the time of interaction between the acoustic 
wave and light. N is called the time-bandwidth product of a particular Bragg 
modulator. [3] 

II. Total Power Time integrating Radiometer 

An AOS can be regarded as a total power, multichannel radiometer with each 
channel modeled as shown in Fig. 3 by a square law detector, with a passband 



o 



o 


filter of bandwidth B, and an Integrator with time constant T. The advantage of 
such a system comes from a processing gain In signal to noise which derives from 
the analog time Integration (cf. Kraus [4], Lucas [5], Kellman et al . £62, and 
a short discussion In Turpin [7]), Specifically, If the input to the square law 
detector Is a Gaussian zero mean signal and Is assumed to be band limited by the 
pre-detection amplifier section, then the filtered voltage, Vjp> resembles a 
randomly modulated carrier wave at the IF frequency. The amplitude distribution 
of the wave form has a Rayleigh distribution given by (see discussion in Kraus 
C43), 


P(v) 



e -(v2/ 2 v2 eff ) 


(U) 


The statistical properties of equation (11) can be found in Whalen [8] and are 

f 

given by 




o “ V 


eft' 


( 12 ) 


that is the mean 
deviation, 0 , of 
viewed as taking 
Integration time 


voltage of equation (11) is equal to the rms standard 
the Rayleigh distribution. The output of :he integrator can 
N * BT independent samples of the diode output and after an 
T the integrator output is given by 


be 


v out * BT “ . 


( 13 ) 


while the rms standard deviation of the output voltage is given by 


o 



<*out * 


/ST 


a 


(14) 


We see that the output of the Integrator, therefore, will be a DC voltage level 
effectively determined by the system noise. If a small signal component Is 
present in addition to the Gaussian noise, the band-limited probability 
distribution In equation (11) evolves to a series of distribution curves called 
Rician density functions (Whalen [8]), with the functions approaching a Gaussian 
distribution for higher input signal to noise ratio. 

The Integrator output with a signal present is given by 

V s - BTS (IS) 

where S is the signal power. Since we are only interested in the signal 
contribution, we subtract off the DC bias of the noise output to derive the 
signal to nc-i M ratio, 

AV out BT S 

SNR 0Ut « - - M SNR in (16) 

°out ^ a 


By integration we obtain an improvement in signal to noise ratio given by /BT, 
which is called the analog processor gain. 

Further sensitivity improvements can be achieved by digital integration 
with extended precision. We digitize the integrator output with an analog to 
digital converter (ADC) after the analog integration period T. With 0 
Integration periods the digital sum is given by the mean value 

“digital 3 J v out 3 J “ BT 

^digital 3 J AV out 3 J s BT ( 17 ) 
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The variance on the digitized mean is given by 


“digital * ^ 0 

The SNR improvement using digital processing Is given by 
SN,! d1g1ta1 ‘ X SNR out ■ 'SBT SNR 1n , 


( 18 ) 


(19) 


and, therefore, we gain a factor /$ improvement in signal to noise ratio over 
the analog process; this gain is only limited by digital precision. 


III. Degradations in Signal to Noise Ratio from Ideal 

In the calculations of analog and digital processing gains, numerous 
additional sources of noise have been neglected which degrades the ideal SNR 
given by equation (19). Among the factors which degrade receiver performance 
and are common to all square-law detector systems are the following: 


a . Recei ver R uctuati ons 


The total system noise is determined primary by the mixer and IF 
pre-amplifier noise figure. Additional stages of gain, especially in a 
room-temperature astronomical heterodyne system which may encompass as much as 
90 dB of amplification, contribute negligibly to the system noise performance. 
This relationship is shown by 


T t Tz Jz 

”sys “ ^l * + + “• 


G z g 2 S 3 


( 20 ) 


where T^ is determined by the mixer noise temperature and conversion losses, T 2 
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Is the noise temperature of the first amplifier with gain given by Gg, and 
further stages of amplification is represented by a noise temperature T 3 , and 
gain G^. Fluctuations In local oscillator power or in mixing efficiencies 
contribute directly to the system temperature. Heterodyne receivers employing 
laser local oscillators and photoconductor mixing elements are particularly 
vulnerable to this type of fluctuation since mixing efficiency and system noise 
depend critically on laser power stability and optical alignment. If the first 
stage of amplification has particularly high gain and Is dominated by Tp 
amplifier stability will also be critical for high sensitivity. 

Receiver fluctuations ulter the sensitivity by adding a post-detection 
noise term to equation (14), and is given by 

« 2 out - BT » 2 RF + < 4|: > 2 (21) 


where we denote a R p aF as the noise contribution from receiver fluctuations, and 
aF is the scale of fluctuations with respect to cr R p. If the receiver 
fluctuations are confined to low frequencies, this is true for amplifiers and 
regulated power supplies, then fluctuation noise can be reduced by synchronous 
detection. In Dicke switching [9], the signal is modulated in a square wave 
fashion against a load, with the modulation frequency chosen to be higher than 
fluctuation frequencies. The signal is recovered by multiplying the detector 
output by a square wave which effectively removes the DC level. However, Dicke 
switching reduces the signal integration time by two, and therefore halves the 
processing gain given by equation (16). The degraded signal to noise ratio 
offered by Dicke switching is given by 


BT 


SNR 


out 


, 7 SNR -in 

/frr+UF-y in 


( 22 ) 


where tJr' is receiver fluctuation noiue reduced by Dicke switching. 


Since almost all astronomical receivers use Dicke switching we will use 


equation (22) for our standard SNR. 


b. Analog to Digital Converter Noise 

The penalty in increasing SNR gain via digital processing is the 
introduction of quantization noise from the finite resolution of an analog to 
digital converter (ADC). In using an ADC, care must be taken to insure that the 
intrinsic noise, aR p, is equal to or greater than the quantization step, i.e. 

a - Ko rf , K < 1 (23) 


If equation (23) is satisfied, then a digitization sample can be modeled to have 
the following statistical properties. An ADC will put out a number proportional 
to the input voltage given by the transfer function iihown in Fig. 4. This 
assumes that the ADC has uniform step sizes and therefore is an ideal step size 
quantifier. The output cm be written as 


V 0 - V i + e (24) 

where e is the error between the ADC outputted code and the true voltage 
value. If equation (23) holds, we can represent the probability distribution of 
e as 


P(e) » (1 /a) for - a/ 2 < e < JZ 


(25) 


Equation (25) implies that the errors are randomly distributed over the 
quantization interval. The mean is given by 
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UNIFORM* STEP SIZE QUANTIZER 
3 UIT MID-RISER 


J»e 


0 


(26) 


+a/2 



(1/A) e 


de ■ 


while the standard deviation is given 


+a/ 2 

«r e 2 - J P(e)(e-u e )' 2 de 

-a/2 


+a/2 


(1/A) 



de 


o + a/2 « ? ? 

- (1/a) e J /3 | » AVI 2 * n c„p /12 (27) 

-a/2 RF 


Equation (26) guarantees that after averaging enough ADC samples we get the true 
average value, i.e. since u e = 0 

7 out “ E i v 0 > 3 E { ^ + e } = E { V i ) + u e - 9 (28) 


It is easy to show that for every additional bit of resolution in an ADC, the 
SNR improves by 6 dB. Since the error distribution is uncorrel ated, the error 
power spectrum will be that of white noise. 

Equation (19) and (22) can now be modified to include the effects of 
quantization noise. 
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WJ 



SNR 


digit 


Bt/J SNR jn 

2 /BT + (aF') 2 + K 2 /12 


(29) 


and the degradation in the SNR Is therefore given by 


a -/l + (aF') 2 /BT - K 2 /12BT (30) 

To insure that receiver fluctuations and quantization noise do not degrade the 
receiver performance, requires that 

BT _ (aF - ) 2 and 

BT _> K 2 /12 (31) 

The proper use or' an ADC requires that K £ 1; if K * 1 (i.e. the quantization 
step size is equal to the noise), the requirement that BT 1/12 in equation 
(31) is easy to meet since BT is usually much greater than 1. Therefore ADC 
quantization noise degrades the SNR negligibly under proper operating 
conditions. 


c. Differential Linearity 

In heterodyne receivers where the aim is ?o recover a small signal out of a 
receiver with large noise, it is important that errors in differential linearity 
of the system be as small as possible.. The major source of differential 
non-1 inearity comes from the use of ADCs. These errors, however, can be 
minimized by proper selection of the converters. In Fig. 5, the widths of three 
different ADCs are plotted for 256 consecutive digital codes. Although the 
manufacturer' s published specifications for all three converters stated ± 1/2 
LSB differential linearity, this performance was satisfied in only one of the 
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three ADCs tested. The results given In Fig. 5 are consistent with the tests of 
the same_ model for other units by the same manufacturer. We can see that the 
worse case Is that of the 12 -bit Analog Devices ADC whose differential error is 
about ± 0.75 tSB, whereas the best case Is the 14-bit Analogic ADC whose 
differential error Is ± 0.1 LSB. In the Hybrid Systems ADC the average error Is 
about j.0.25 LSB but with excursions of + 1, - 1/2 LSB in this voltage range. 

Differential non-linearities affect our measurements by altering the 
incremental linearity or the digital transfer function from an ideal one. The 
transfer function of a particularly bad ADC is illustrated in Fig. 6 in contrast 
to an ideal incremental transfer function shown in Pig. 4. The bit sizes, a., 
in Fig. 6 are all different as compared to the uniform step sizes given by an 
ideal ADC in Fig. 4. In an observation, the actual transfer function is 
convolved with a Gaussian of width which is represented in Fig. 6 by the 
smooth curve; however, the convolved digital transfer function will depart from 
the ideal one given by the unit slope dotted line. 

In detecting a small signal by means of subtracting a large DC offset 
representing the receiver noise, the absolute error in the digital transfer 
function is not the determining factor which gives the error in the measurement, 
rather it is the departure from unity of the slope of the transfer function at 
the reference operating point, i.e. since 


A S * »(V) A V in 


(32) 


where m is slope at the reference voltage, V, and aV is the signal difference. 
The error in the measurement can then be given by 


e(V) - A rn .(V) aV 


(33) 
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Incremental Transfer Function 


FIGURE 6. Non-ideal ADC incremental transfer function. 

Solid line represents convolved non-linear 
transfer function compared to ideal ( dotted 
1 





where E (V) is the error function of the reference voltage, ^m(V) is the 
deviation from unity of the slope of the digital transfer function, i.e. 
m(V) - 1 + A m(V). 

Since the error Is a function of V, the errors of independent measurements 
of aS can no longer be uncorrelated. This Is due to the fact that in deriving 
the smoothed digital transfer function, we convolved the incremental transfer 
function with a Gaussian of width a, which defined the deviation we can expect 
from V. Therefore Am(V) will be, on the average, the slope error we can expect 
at the reference voltage. We then see that C (V) will be the limiting 
sensitivity since we cannot integrate to a value below c(V), defined by the 
differential non-linearity error. In the previous sections we saw that with N 
observations of a signal, the signal increased as N, i.e. S * NaS, while the 
uncorrelated noise increase as /fT, i.e. a * s * No^> therefore * A • 
and the SNR 0Ut * s/cr s « /h SNR^ n . However, with correlated noise, the noise or 
error new increases as N, a * Ne(v) and the output signal to noise SNR 0Ut * 

SNR ln =* s/e(v) which is limited by the differential non-linearity error of the 
transfer function. 

It is obvious that differential non-linearities in any of the amplifiers 
which are in line before the detection stage will have the same consequences. 

If an amplifier is operated near the saturation range for a particular frequency 
component in its bandpass, the differential non-linearity errors may be quite 
large. 


d. Calibration 

Calibration of multichannel data serves two purposes, the first is to 
remove channel by channel variations of gain in the system, and in so doing 
accomplishes the second which is to establish a physical scale for the 
measurements. Calibrating an observation can degrade the SNR of the 
observation. Although this degradation is not an instrumental effect, a 
detailed examination of the calibration process can show us how to minimize this 
type of degradation. 


IS 


u 


Schematically, the output, Sj, of each filter channel is given by 


s < - G 1 < T sys + C 1 s 1> 


(34) 


where G^ is the channel by channel gain, 1 - 1 to N where N is the mber of 
channels In a fllterbank, T sys Is the receiver system temperature, s i is the 
signal Input to the receiver, and represents the frequency response of the 
receiver front-end. The distinction between and ei Is important. Although 
the Gj's are different channel to channel, they are fairly close in value within 
the entire IF bandpass; l.e. we choose the RF power amplifiers to have 
sufficient bandwidth. The front-end response, , however, may have large 
variations across the IF bandpass. This is especially true in the lOum 
heterodyne receiver (Mumma et al . [10]), where falls in sensitivity by a 
factor of 3 from 0 to 1.6 GHz in the IF. In addition, depends strongly on L0 
power, optical stability, and mixing efficiency. It is Important to note that 
even in the absence of a signal, , there will st 11 be an output * Sa T sys • 
Therefore, the front-end gain and time stability can only be sampled if a 
sufficiently strong signal is present, i.e. e-- comes only into play when an 
input signal, s^ , is not zero. 


The calibration scheme works the following way. A measurement is taken of 
a known temperature blackbody, T BB , usually chopped with reference to cold sky. 
Although the sky temperature is 270°K, observations are undertaken where the sky 
optical depth is relatively small, i.e. where the sky is not opaque. Therefore, 


it is a good approximation , for our case, to t«ke T 


sky 


0. We get 


s i * ^ ( T sys + e i T BB ) 

*1 * G 1 ^sys + c 1^sky) “ G 1^sys 


We take the ratio (S^R^J/Rj channel by channel 
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Sj-fy G i( T sys * M T B B ) " G 1 Ts y s 

m SI 

R 1 S 1 T sys 

e i t BB 

■ 

^sys 


( 36 ) 


where BB^ is the blackbody measurement, and T sys Is generally not known £ 
priori . BB^ is actually a good measurement of 1/ Ci which Fig. 7 shows very 
plainly. A similar ratio Is taken with a signal at the receiver input given by 
aTj, The measurement, » is given by 


aT 4 




sys 


(37) 


Now we can solve for aT^ » using equation (36) since we can scale with the 
blackbody temperature Tgg. (T BB is not measured exactly through the same optical 
path as the signal, l.e. we will ignore telescope or transfer optic efficiencies 
here.) 


M 1 

* T 1 * T BB ( 38 > 

BB i 

and we have eliminated all gains from the measurement. 

Assuming BB^ , and T BB are noiseless measurements, the only contribution in 
uncertainty in measuring A T^ comes from the uncertainty in measuring M-j . We can 
calculate the SNR of a measurement of the type M * (S-R)/R. The noise 
contribution to M is given by 
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Fig. 7 . Measurement of system degradation A as a function of IF 
frequency’ (0 to 1600 MHz). The three-fold Increase In A at high IF 
frequencies Is caused by. the roll-off. of. the mlxer-preampllflai 
response.' • ’ • 




Using a s ■ or we get 



(40) 


The signal to noise ratio for M Is therefore given by 



C M 



(41) 


where we acted that (S-R)/ a is the SNR of the signal itself given by equations 
(19) and (£0), the signal to noise ratio has been degraded by 


/ 2 

& » /l + (S/R) 


(42) 

Let S = aS+R, where aS Is the signal response over the reference. We get from 
equation (40) that 


When R»aS, then S - Ji. % e.g. aS» 0 In a balanced receiver mode; the /l SNR 
degradation results simply from taking the ratio. For S/R > 1, then a * S/R 
from equation (40). 

It is now easy to see why a Dicke-switched balanced receiver has the 
highest sensitivity. A balanced receiver operates at the point S ^ R, and 
aS«R. In absorption observations of atmospheric constituents against the sun 
or the moon, aS can become a sizable factor with respect to ft, especially if the 
reference is against the sky, Tc minimize the degradation in SNR of solar 
absorption spectra, the observations should be chopped against a reference 
temperature which is greater than, or equal to the sun's, i.e. R £aS, which means 
S/R * 1. 

IV. Photodiodes as Square-Law Detectors 

In an AOS we must Include noise terms which are unique to using 
photodetectors as total power Integrators. These terms are in addition to the 
sensitivity consideration given in Section III which are common to all total 
power radiometers. Within the AOS, the photodiodes in a linear array acts as 
the basic power detectors - they serve the functions of bandpass filters, square 
law detectors, and power integrators. Since the intensity of diffracted light 
is proportional to the spectral power of the IF (given by equation 6) at a given 
frequency, a photodiode, at a particular location, integrates the light and 
gives an output which is the integrated power at that frequency component. Each 
output of a pixel in an AOS linear array, therefore, acts as the single channel 
as shown in the schematic of Fig. 3. 

We will discuss here, primarily, a set of commercially available linear 
arrays of monolithic silicon integrated photodiodes with 1024 elements. Each 
photodiode is on a 15um center with dimensions of 7um width by 16um (or 300 urn) 
heights. The diode has an associated storage capacitor on which to integrate 
photocurrent and a multiplex switch for periodic readout via an integrated shift 
register scanning circuit. Some properties of a typical diode array (Reticon 
1024H) are shown in Fig. 8 [II^ Fig. 8(a) is the detector spectral response, a 
typical quantum efficiency at the HeNe laser frequency of 632.8 ran is n = Q.7, a 
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CW GaAs diode laser operates at \ » 850 nm, with n 0.7. The geometry and 

aperture response of a pixel Is shown in Fig. 8(b). The equivalent circuit of 
the array is given In Fig. 8(c); notice a row of dummy diodes In parallel to the 
light sensitive ones. 

a) Detector Responsivity and Storage Capacity 
The responsity of a photodetector is defined by 


9nU) 
Rx a 

Ahc/x 


2 

(in units of ampere/ watt/m or 
coul ombs/joul e/nr) 


(44) 


This number is dependent on the wavelength of light through the quantum 
efficiency n(x)» and the photon energy hc/x; q is the electronic charge and A is 
the photoelement exposure area. Typical values for broadband responsivity (over 
the dotted light source given in Fig. 8(a)) for the Reticon 1024H is R*0.8 
pA/^watt/cm^. In order to convert equation (44) to the number of photoelectrons 
produced by an exposure to light we have 


n e (T) » (R x /q)(hc/ x )(* x AT) = R^T 


(45) 


where * is the flux given by number of photons/sec, and T is the exposure time 
for the photoelement, and R Qx is defined by (45). 

Each photoelement has a finite charge capacity which is given by its 
associated on-chip capacitor. For the Reticon 1Q24H tne storage capacity is q = 
2.2 pcoul which corresponds to a electron capacity of N sat 3 1.37 x 10^ 
electrons. The maximum integrating time for an array, in the absence of 
detector noise, is then given by 
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T < N sat /( n *) 


(46) 


b) Dark Current and Leakage Current Noise 

A photodetector, even in the absence of light, will generate current due to 
the diffusion of charges generated by thermal excitation. This leakage current 
will charge up a storage capacitor and, since the generation process is random, 
contributes temporal noise associated with the sensing of the charge or voltage 
on the capacitor. The process obeys Poisson statistics and can be described as 
shot noise. The error given by the process (in units of electrons) is given by 


Q n 2 (leakage) * qijJ 


(47) 


where i^ is the dark or leakage current and T is the sample time interval. Even 
in the absence of light, therefore, the leakage current will limit the maximum 
integration time 


^max 1 ^sat^ L 


(48) 


For the Reticon 1Q24H, ijj»1.0 pA at room temperature. The leakage current is a 
very strong function of temperature and will double in value for every 7°C. The 
dark current contribution will saturate the detector in T mav =,2.2 sec at room 
temperature (25°C), 

c) Charge Reset Noise 

The noise contribution from the charging of capacitors is important since 
in a photodetector array this mechanism is used repeatedly to convert a signal 
from charge to voltage and vica versa (see Borsuk [12]). Within a detector a 
capacitor whose initial charge is zero is to be charged to a voltage V Q , through 
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a resistor R by closing a noiseless switch. The voltage across the capacitor is 
then given by 


V(t) * v 0 (1 - e" t/RC ) 


(49) 


or in terms or' the charge 


l)(t) - CV 0 (1 . e‘ t/RC ) 


(50) 


The variance on the mean value or' q(t) has been shown by Barbe [13] to be given 
by equation (51) where T is the temperature and k is Boltzmann's constant; 


Q(t) 2 - kTC (1 - e ' 2t/RC ) 


(51) 


The signal detection process involves two steps. The first is the fast charging 
up of the capacitor by the action of a low resistance switch. The second comes 
from the discharging of the capacitance by photogenerated charges when the 
switch is turned off. As noted before, thermally generated leakage current also 
contribute to this process and adds shot noise. The RC/2 time constant 
associated with equation (51) is different for the two processes. In the 
charging mode a typical MOS switch, which has a nominal value of R * io^q 
charging up to 2.2 pcoul , has a time constant RC/2 * (10 4 ax0.2pF)/2 = 1 nsec. 
When the switch is turned off, then R * 1Q 12 q, and RC/2 « 0.1 sec. 

Consequently, after t * SRC in the charge mode the mean deviation is given by 


% - (kTC) 


(52) 


i/2 


If we sample the charge immediately after the low resistance switch has been 
turned off, we can then get a measurement of Q n without any additional noise 
since the time constant Is so long. This noise contribution can be taken out by 
subtracting it from the measurement after an exposure time t. 

This type of signal processing ■‘s commonly called correlated double sampling 
(CDS), and is used to subtract out the charge reset noise. We can calculate the 
noise equivalent electron number given by equation (52) since fl e * Q n /q 3 
✓kTC/q. In terms of C given in pF and T 3 300°K we get 

R e - 400 /T (53) 

The detector storage capacitance on the Reticon 1024H is given by C * ^ sa tV v 0 
where Vo is the charging voltage. q sat * 2.2 pCoul and V Q 3 12V, therefore C * 
2.2 x 10"^/12 * 0.18 pF. The noise equivalent electrons is then given by 
equation (53) R e * 400 /Q.18 <* 171 electrons. However, each detector is then 
read out of a common video line which has a 35 pF capacitance. This readout 
noise is, N reac j 0ut * 400 ^ a 2336 electrons. Therefore, the mean noise is 
dominated by video readout noise. It must be noted that the charge reset noi.'je 
is a post -detect ion factor which will occur at the same level after the sampling 
of each detector element, and is therefore independent of the exposure time. 

d. Signal and Background Shot Noise 

Whenever discrete particles arrive at random times, there will be 
fluctuations in the rate of arrival.. The detection of light involves the 
measurement of the arrival of photons,. and therefore the fluctuations in the 
rate of arrival constitutes an additional noise source which is called shot 
noise. This type of process can be analysed in a straight forward way (see 
Oliver [14]) and we will briefly- derive the result. We assume that a photon is 
equally likely to arrive at any time but that it arrives at an average rate, r. 
This means that if we measure n photons arriving within a time t, the limit as t 
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* » of n/t * r. In a time interval at « 1/r, the probability of detecting a 
photon is just rat and the probability of not detecting a photon is (1-rat). 
Therefore the probability of no arrivals during t is given by 


p (t) - lira (l-r A t) t/At 
a at+o 


(54) 


since there are t/at independent intervals during which no arrival must occur, 
and taking the limit, equation (54) becomes 


P 0 (t) - e‘ rt (55) 

The probability of exactly one arrival in time t is the probability of receiving 
a particle in the Interval T to t+dr and none before or after, summed over all 
f • We then have 


z 

, l( t ) * J P 0 ( t ) C rdT 3 P 0 ( t - T ) 


3 rt e 


-rt 


(56) 


For the arrival of two photons, the probability is given by one arrival before 
x, one arrival in the interval r+d-r, and none after, i.e. 



. (rt) 2 e -rt 

71 


(57) 
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And In genera] tor n arrivals we get 


Pn«) 


(rt) n -rt 


n! 


The mean value of n f is given by 


(58) 


n » z n p (t) » rt (59) 

» 

from evaluating equation (58) and (59), as It should be from the definition of 
rt. We can obtain a recursive definition for p (t) from equation (59) and get 

(n) n _ n 

P(n) - e n * „ p(n-l) (60) 

nl n 

The deviation from the mean is derived from 

(an) 2 = z P(n) (n-n) 2 

* n 2 - (n) 2 (61) 

and we can calculate 


n 2 * z n 2 p(n) ■ fi + (n) 2 . 


( 62 ) 


Therefore we get the simple result that 
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0 



(an) 2 » fi or an a A (63) 

In an AOS, the diffracted light level detected by the photodetector array will 
be a large constant, due to the receiver noise, plus a small contribution due to 
■ « the signal, in addition a contribution due to scattered ambient light from the 
laser can be present. If we represent the photon flux from the receiver noise 
plus signal as * s and the background term as then the detector output Is 
given by equation (45), and the noise contribution will be given by 


(sn) 2 - R 0x (« s * « b )T (64) 


Careful optical arrangement can minimize * B , but the contribution from the 
signal term will always be present. 

If we operate near saturation n g « l x 10 7 electrons, then the dynamic range is 

limited by the signal shot noise as given by nj/lT « 3 x 10 3 or 30 dB. We see 
3 e e 

that jn ■ 3 X 10 electrons, and this value is comparable to the charge reset 
noise of * 2 x 10 3 electrons derived in the previous section. 

c) User Stability and Mechanical Alignment 

Two potential sources of receiver fluctuation inherent in an AOS is the 
stability and added noise of the laser light source. In a bulk Bragg cell 
system a HeHe laser Is used as the coherent light source. Since the optical 
components are made up of discrete parts, overall mechanical stability of the 
system must also be considered* In an integrated optics AOS, the laser light 
source is a solid state GaAlAs diode laser which is precisely butt-coupled to 
one end of Li NbQ 3 substrate on which lenses and transducers are fabricated. 
Mechanical alignment and jitter have been minimized by the compact and 
structurally integral package. However diode lasers may be susceptible to ion 
plasma noise. 

The fluctuations in power of the Spectra Physics Model 120 HeNe laser, 
which we use in our bulk Bragg cell AOS, is specified to be <0.5% from 1-100 
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kHz, and amplitude ripple Is less than <0.2% up to 120 Hz. The primary cause of 
amplitude noise in the HeNe laser is intermodulation of the optical frequencies, 
primarily five optical resonant frequencies spaced about 300 MHz apart within 
the 1500 MHz half-width temperature-broadened neon gain profile. When the 
resonant frequencies are combined in a non-linear detector, as In a photodiode, 
the output contains all possible difference frequencies between the optical 
frequencies present. These intermodulation products appear noise-like since this 
frequency components vary rapidly with small cavity changes. The 
intermodulation products limit the minimum noise level in the 1 KHz to 100 KHz 
range of the HeNe laser. 

A photodiode in addition to acting as a square-law detector Is also an 
integrating element, and therefore acts as a low pass filter. For an 
integration time in a photodetector of *10 msec, the detector filters out 
fluctuations above 50 Hz (the Nyquist frequency). Fluctuation of the 5 mW laser 
at this frequency has been specified at 0.2%. Noise sources at these 
frequencies can be treated as gain fluctuations as described previously at in 
section Ilia, and can be minimized by Dicke-switching. 

In a bulk Bragg cell AOS system, the primary consequence of mechanical 
jitter Is the frequency stability of the spectrometer caused by errors in the 
optical pointing. Proper mechanical design and rigid mounting of optical parts 
can minimize this problem. Measurements of frequency error in a prototype AOS 
system has shown that optical alignment can be attained to within a detector 
pixel, about 10 u m. Fig. 9 shows the frequency stability of a single tone 
observed on an AOS over a timescale of over 3 hours. The ordinate scale 
represents channel number or position of a single pixel. 

A potential source of error in using an AOS for detecting signals which are 
as narrow as the frequency resolution of an AOS is an effect called "scalloping 
loss." This effect results from the fact that output plane of the Bragg cell is 
sampled only at discrete positions with an aperture resolution which does not 
totally overlap neighboring detectors (see Fig. 8), i.e. if a narrowband 
frequency component happens to fall between two detectors, the resulting 
amplitude output will be divided between the two pixel elements. In this case a 
measurement of the strength of the signal will be seriously underestimated. 
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This problem Is entirely analogous to the use of digital FFT's to estimate 
spectral components when the number of points In the FFT is not large enough to 
sample the narrowest frequency component. This effect is usually not a serious 
one in astronomical applications, since spectral features are often much wider 
than a single resolution element. 

1 l 

V. Summary of Noise Terms 

We see in section II that an Improvement In SNR from the original input SNR 
can be obtained by both analog and digital processing. This is an essential 
feature since the input SNR can be as low as 10“ 3 or 10" 4 with submillimeter 
heterodyne receivers now being used or under development. However, degradations 
to this pro cessing gain comes from several sources; some outlined in section III 
are common to all total power radiometers, and others as outlined in section IV, 
are uniquely attributed to the AOS which depends on photodetectors and a 
coherent light source as the detection system. We rewrite equation (22) with 
the degradation term denoted by a. 


✓5bT SNR, n 

SNR out - IL 

2 6 (65) 


where BT Is the time bandwidth factor, and d comes from the number of digital 
samples. T is limited in an AOS by several factors, primarily by the saturation 
capacity of a photoelement. J, however, is only limited by the digital 
precision in which the digital processing takes place. The degradation can be 
written as 


6 


2 


1 + 


N 

z 



( 66 ) 


where we denoted the extra noise contributions by 


2 


which we sunmarize below. 
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1) Gain Instability as described In section Ilia and which may Include i-ser 
and mechanical stability described In section IVe, This term Is given by 


■ (iF ') 2 /BT 


(67) 


where aF' percentage of <j RF gain noise present even after Dlcke-swltch 
filtering, 

2) ADC (juanti zation Noise 


8 2 2 - K 2 /12BT 


( 68 ) 


where K<1 and K Is the scale factor of the quantization step with respect to the 
RF noise, l.e. a * Ko RF . This term can be Ignored in a properly operating 
system. For example K » 1 using the worse case, B * 1 MHz, T » 1 x icf 3 sec, 


then BT » 10 , which means 


« 2 2 * 8 X 10~ 5 


3) Diode Leakage Current 


«3 * ’*L^ <r RF B ) 


( 69 ) 


4) Diode Charge Reset Noi 


KTC 


16 x IQ 4 C 


q2ff2 RF Bt 


or nr-Bt 


RF 


( 70 ) 
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where C is measured In pF, and where t Is now the Integration time, and T is the 
temperature. 

5) Signal and Background Shot Noise 






2 R 
a rf B 


(71) 


The above five terms are uncorrelated noise terms which do not alter the fact 
that the processing gain increases as A)BT given by equation (65). However, we 
also pointed out in section 1 1 Id, that the errors associated with differential 
non-linearities , are correlated with one another, and therefore cannot be 
processed out. The processing gain may be limited by 


SNR 


out 



(72) 


where S/c(v) is the SNR determined by the differential non-1 inearity at a 
particular operating point. 

Finally, we have a degradation in SNR due to the arithmetic processing of 
the observed data In order to cancel channel to channel gain variations. 

Although this effect is not instrumental in nature its analysis enables us to 
minimize this factor by reconsidering our observing scheme. This degradation is 
discussed in section I lie, and is an additional term which must degrade equation 
(61). 
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Process 


X 


SNR 


SNR 


out 


« 

process 


where 


4 process 


1 + (S/R) 2 


(73) 


( 74 ) 


where S is a measurement of a signal and R is a measurement of a reference load 
In a Dicke switch system. 
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